Ac-cotunneling through an interacting quantum dot in a magnetic field 
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We analyze inelastic cotunneling through an interacting quantum dot subject to an ambient 
magnetic field in the weak tunneling regime under a non-adiabatic time-dependent bias- voltage. Our 
results clearly exhibit photon-assisted satellites and an overall suppression of differential conductance 
with increasing driving amplitude, which is consistent with experiments. We also predict a zero- 
anomaly in differential conductance under an appropriate driving frequency. 
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Recently, cotunneling^ through discrete levels, i.e., 
a quantum dot (QD), has attracted much attention 
since it determines the intrinsic limitation of accuracy 
of single-electron transistors due to leakage, and since 
it also involves correlation effects, such as the Kondo 
effect^ It has also been reported experimentally and 
theoretically^^ 7 - that external microwave irradiation 
can induce the occurrence of Kondo satellites and an 
overall suppression of the Kondo peak. However, there 
are few studies so far concerning time-dependent second- 
order cotunneling in the weak tunneling regime at tem- 
peratures above the Kondo temperature. About ten 
years ago, Flensberg presented an analysis for coher- 
ent photon-assisted cotunneling in a double-junction 
Coulomb blockade device in the adiabatic limit.— In this 
letter, we will further study the cotunneling in an inter- 
acting QD when an ac bias-voltage is applied between 
two electrodes in the non-adiabatic regime. 

We employ the s-d exchange Hamiltonian to model in- 
elastic cotunneling through a QD in an ambient magnetic 
field, B, in the weak-coupling regime: 
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where c^ k(J (c^ko-) is the creation (annihilation) opera- 
tor for electrons with momentum k, spin-er in lead 77 
(= L,R). The energies e v k(t) = e° )k + eV v (t) include 
a rigid shift of the Fermi energy of the electrons in 
the leads due to the applied time-dependent bias-voltage 
V v {t) = V2 + v v cos{VLt) with VS (v n ) being the ampli- 
tude of the dc(ac) part of the bias-voltage. Here, we as- 
sume that the Fermi energies of two leads are zero at 
equilibrium and VjP = —V§ = eVo/2. A = g e HBB 



is the static magnetic-field _B-induced Zeeman energy. 
S = (S x , S y , S z ) are Pauli spin operators of electrons 
in the QD [S^ = S x ± iS y ], and J is the exchange cou- 
pling constant. -f/dir is the potential scattering term with 
2Jd = J. As in our previous paper p we can rewrite the 
tunneling term, Hi in Eq. (fT]), as a sum of three products 
of two variables: 



Hi = Q Z S Z + Q+S- + Q-S+ + Q\ 
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with the generalized coordinates Q z ^ of reservoir vari- 
ables as 
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and Q 1 — -ffdn- In the following, we will use units where 
H = ks = e = 1. 

As in our previous studies of inelastic cotunnel- 
ing through an interacting QD in the weak tunneling 
limit, we employ a generic quantum Langevin equation 
approac h 10 ' 11 ! 12 to establish a set of quantum Bloch equa- 
tions for the description of the dynamics of a single spin 
[modeled by Eq. ([I])] explicitly in terms of the response 
and correlation functions of free reservoir variables. This 
procedure provides explicit analytical expressions for 
the nonequilibrium magnetization and cotunneling cur- 
rent for arbitrary dc bias-voltage and temperature^^ 
Here, we generalize our previous derivations to the time- 
dependent case in the non-adiabatic and high-frequency 
regime. 

In the derivation, we proceed with the Heisenberg 
equation of motion for the Pauli spin operators and the 
lead operators, and then formally integrate these equa- 
tions from the initial time to t exactly to all orders of 
J. Next, under the assumption that the time scale of 
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decay processes is much slower than that of free evolu- 
tion, we replace the time-dependent operators involved 
in the integrals of these EOM's approximately in terms 
of their free evolutions. Thirdly, these EOM's are ex- 
panded in powers of J up to second order, resulting in 
non-Markovian dynamic equations for the time evolution 
of the QD spin variables in a compact form: 

S z (t) = -2 f dt' {e~ iA ° T + e iA ° T )C{t,t')S z {t') 
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with t = t — t! . The correlation function, C(t,t'), and 
the response function, R(t,t'), of free reservoir variables 
(tagged by subscript "o") are defined as: 

C{R){t,t>) = if(r)<[Q±(t),Qj(0]+(-)>. (8) 

Special attention must be paid to the free reservoir 
variables due to the time-dependent energies e^k(i) : 
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Therefore, the kernels C(R)(t,t') become double-time- 
dependent functions due to the lack of time-translation- 
invariance stemming from the ac-bias: 

C(R)(t,t') = \e{T){[Qt{t),Q-{t')] ± ) 
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with g = J 2 Pq, V ac = vl — vr, and the Fermi-distribution 

function is f n (e) = [l + e^ £ ~ M "-'/ T ] (T is the tempera- 
ture). Here we assume the two electrodes to be Markov- 
type reservoirs with a constant density of states pq. The 



kernels reduce exactly to our previous results, Eq. (B8) 
in Ref. [9j, if there is no ac-bias or with the same ac am- 
plitude in the left and right leads. 

In the presence of a periodic ac-bias, the spin variables 
S z ^(t) naturally depend periodically on t with a period 
T &c = 2tt/Q, As a result, the full solutions of Eqs. (J6j) 
and |7]) can be formally written as a superposition of all 
harmonics 
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Employing this expansion in the dynamic equations §6§ 
and ([7]) , an infinite set of linear equations results in which 
the o A ' are coupled with each other via the kernels. To 
obtain a solution for the spin variables, one has to termi- 
nate this infinite chain at a chosen order and then solve 
the resulting equations in a recursive way. However, in 
the limit of high frequencies, f2 ^> T (tunneling rate) 
and T (temperature), of interest in this letter, the ac- 
bias oscillates so fast that an electron experiences many 
cycles of the ac-bias during its presence inside the dot, 
and thus can not sense the details of the dynamics within 
one period T ac . In this non-adiabatic limit, one can ap- 
proximately replace the kernels by a time-average with 
respect to the center-of-mass of time f = t + t'i 14 ' 15 
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By the same token, one can retain only the station- 
ary part of the spin variables and neglect the rapidly 
oscillatory parts, leading to detailed balance equations 
in a Markov approximation by making the replacement 
r*^ dr => dr in Eqs. © and Q: 
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Here, C(i?)L 0) H are the Fourier transforms of the time- 
aver ag ed kernels C{R)^(t); 
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with (f(x) = x coth(x/2) and J n (x) is the Bessel function 
of order n. To derive these equations (fT3)) - (fTrJ|) , we use 
the relation 



E Jn ^ e 
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The solution of Eq. ((T3]) yields the nonequilibrium 
magnetization of the QD subject to an ac-bias voltage 
as 
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This formula is our central result, which can be re- 
garded as a direct generalization of the dc nonequilib- 
rium magnetization ^ 16 ' 17 of a QD under a non-adiabatic 
high-frequency field. Obviously, it reduces exactly to pre- 
vious results in absence of ac-bias, V ac = 0i£ i 16 i 17 As an 
illustration, we exhibit in Fig. 1(a) the dependence of 
the magnetization, Sf \, on dc bias- voltage for the driv- 
ing frequency f2/Ao = 0.5. It should be noted that Sfo 
exhibits different behaviors with increasing ac-amplitudc 
14c- For small dc bias voltage, the QD spin is fully po- 
larized due to the nonzero external magnetic field, and 
it is gradually quenched with increasing dc bias voltage. 
Application of an ac bias tends to quench the spin po- 
larization more rapidly. This tendency suggests that the 
ac bias plays a role in dephasing the electronic tunneling 
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which is a generalization of Tien-Gordon-type formula in 
cotunneling current.— It should be noted that our ex- 
pression for photon-assisted cotunneling current is valid 
in the high frequency limit, whereas Flensberg derived a 
perturbative current under pump conditions in the low- 

frpmipnpv limit § 
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Figure 1: Nonequilibrium magnetization, Sf s (a); and the 
differential conductance, dl/dVo (in units of 2ge 2 /h) (b); as 
functions of dc bias- voltage, Vo, for various amplitudes of the 
ac-bias with a fixed driving frequency fl/Ao — 0-5- The tem- 
perature we use in the calculation is T/Ao = 0.01. 

Proceeding to the calculation of tunneling current, the 
current operator through the QD is defined as the time 
rate of change of charge density N v = J^k.a c lka c r,k<? m 
lead n: I n (t) = N n . From linear-response theory we have 

I(t) = (I v (t)) = -i [ dt'([J v (t), iTi(t')]->- (19) 



Because the dc component of current is easily measurable 
experimentally, we compute the time-averaged current 

1 = %r /o dtI ^)- Perfor 

ming the same high-frequency 
and Markov approximations as above, we obtain the 
time-averaged currents: 
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Figure 2: Calculated differential conductance vs. dc-bias. (a) 
results for a relatively high driving frequency fi/Ao = 1.5 
for various ac-amplitudes, as in Fig. 1; (b) the ac- frequency 
dependence of dl /dVo for a fixed driving amplitude Vkc/Ao = 
1.0. Other parameters are the same as in Fig. 1. 

We plot the dc bias- voltage-dependent differential con- 
ductance, dl/dVo, in Figs. 1(b) and 2. Obviously, the 
differential conductance shows some satellites at Vo = 
±(Ao — nfl) superimposed on the characteristic jump at 
Vo = ±Ao in the presence of an ac-bias. These satel- 
lites arise physically from photon-assisted spin-flip co- 
tunneling, i.e., albeit |Vb| < Ao, the spin-flip cotunnel- 
ing process can still become energetically activated by 
an electron absorbing photon quanta to compensate for 
the energy difference. Moreover, an overall suppression 
is observed with increasing ac-amplitude, which is qual- 
itatively consistent with the experimental results^ More 
interestingly, we find that dl/dVo exhibits a transition 
from peak-splitting to zero-bias- anomaly if the driving 
frequency is higher than Ao. This behavior can be as- 
cribed to the fact that a driving field with an appro- 
priately high frequency can spur spin-flip cotunneling, 
notwithstanding \Vq\ < (fi — Ao); in contrast, when the 
dc-bias increases to Ao > |Vo| > (fi — Ao), spin-flip 
events become inactive instead. The dl/dVo curve re- 
covers peak-splitting behavior if > 2Ao, as shown in 
Fig. 2(b). 

In summary, we have generalized the generic Langevin 
equation approach to study photon-assisted cotunneling 
through an interacting QD in the non-adiabatic and high 
frequency regime, deriving explicit analytic expressions 
for the dc components of nonequilibrium magnetization 
and current with a generalized Tien-Gordon-type form. 
Our results show that applying an ac-bias is an important 
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method for tuning the I-V characteristics. Considering 
experiments^ in which a static magnetic field B = 11 
T is applied, A ~ 0.1 meV, and the ac-frequency is 
Q ~ 12 - 50 GHz with ac-amplitude V ac < 0.15 mV, all 
the parameters are easily accessible experimentally and 
they satisfy the non-adiabatic condition Q, ^> T (~ 30 



M eV). 
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